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Abstract 

The problem of stochastis in multiple electron-photon head-on collisions has been 
considered in this paper. The kinetic equations for the distributions over the elec- 
tron energy and collisions number along with the equations for these distributions 
moments have been obtained. The equations for the first moments have been solved 
by the iteration method. It has been shown that the variance of the energy distribu- 
tion as a function of the mean number of collisions has a maximum at some value of 
n. It is seen from this analysis that multiple scattering of electrons influences on the 
spectra both electrons and photons even for the photon target of small thickness. 
The data of approximate analytical calculations agree with the results of the Monte 
Carlo simulation. 
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1 Introduction 



The problem of laser light Compton scattering on high-energy electrons are 
now considered in the projects relating to the creation of 7 — 7 colliders, 
laser-synchrotron sources, laser cooling, diagnostic of sub-picosecond electron 
bunches and others [1]. It is supposed that intensity of laser flash in these 
problems is so high that an electron can undergo several successive collisions 
passing through a photon bunch [2-4] . The distributions over the electrons en- 
ergy and collisions number, the moments of these distributions, and the spec- 
tra of scattered photons have been studied in this paper using corresponding 
kinetic equations and by statistical simulation methods. It is seen from this 
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analysis that multiple scattering of electrons influences on the spectra both 
electrons and photons even for the photon target of small thickness. 



2 Kinetic equations 



Penetration of an electron through a photon target is a stochastic process, 
where both the number of collisions of each electron with photons and the 
energy loss in individual collision are random. 

Let us consider an electron with energy e traveling through a bunch of pho- 
tons with energy uj . The typical value of the scattering electron angle in the 
Compton back-scattering process is determined by the laser photon energy u 
and doesn't depend on the electron energy e : 



mc 



mc 2 being the rest energy of electron. We shall consider the head-on collisions 
of laser photons with energy ujq ~ 1 eV and electrons with energy e ~ 10 
GeV. In this case the electron deflection angle is much less than characteristic 
radiation angle mc 2 /eo and one can neglect the electron angular deflection 
(the straight-ahead approximation). In this approximation the probability to 
undergo n collision along the pass /, P(n\eo,l), obeys the adjoint balance 
equation (the Kolmogorov- Chapman equation) [5,6]: 



P(n|e ,Z) = ( 1 - sE(e ))P(n\e ,l- s) 

Wmax „ , n 

+ s£( £o ) / -^?-P(n-l\e -u,l)du, (1) 



S(eo) and £(oj;£o) being the total and differential macroscopic cross-sections 
of the Compton scattering, s is a small part of I, and uj is the energy of 
scattered photon (0 < uj < w m ax), 

x 

Wmax = — ; 

1 + X 



is the maximum value of the scattered photon energy, 



X = 



(mc 2 ) 2 ' 
E(e ) = 2n L a(e ) , 



2 



w x da(u\ eg) 
E(u;;£:o) = 2n L 3 



cr and — are the total and differential cross-sections, is the concentration 
of laser photons in a bunch. 

Note that S(eo) is the mean number of collisions of an electron per unit path 
length and Tj{uo\sq) is the mean number of an electron collisions with energy 
loss in unit interval about uj per unit path length. 

The first term in the right side of Eq. (1) corresponds to electrons which 
pass the path s without collisions and 1 — sS(eo) is corresponding probability. 
These electrons have to undergo n collisions along the rest path I — s. The 
second term corresponds to the electrons which undergo the first scattering 
passing the path s and sS(e ) is corresponding probability. These electrons 
have to undergo n — 1 collisions after that but the energy of electron after first 
scattering equals sq — uj, where uj is random energy of the scattered photon 
and £(a>; £o)/£(£o) is the probability density function of uj. In the limit s — > 
Eq. (1) gives the integro-differential equation for P(n\e ,l) [7]: 



d 

—P(n\e ,l) + E(e )P(n\e ,l) 



dl 



^max 

-J Z(u;e )P(n- l\e -w,l)dw = (2) 



with boundary condition 

P(n|e ,Z)|j=o = $no ■ 



In a similar way one can obtain the kinetic equation for the probability den- 
sity function P(e\e ,l) describing the energy distribution of electrons after 
travelling the path I: 



d_ 

dl 



*max 

P(e\e ,l) + S(e o )P(e|eo,0 - / e Q )P{e\e Q - w, l)dw = (3) 



with boundary condition 

P(e\e ,l)\i= = 5(e - e ), 



3 



5(e — e) being the Dirac ^-function. 

Eqs. (2), (3) can be transformed into the equations for the moments of distri- 
butions P(n\e , 0> P( £ \ £ o,l)' 



oo 



(e ,l) = ^n k P(n|e o ,0> 

n=0 
£(> 

e*(e ,l) = J e k P{e\e ,l)de. 



The equation for n and e k has a form 



d_ 



dl 



Wraax 

e Q , I) + Y,(e )n(e , 1) - J S(cj; e )n(e - u, l)dw = E(e ) , (4) 



—e k (e ,l) + Z(e )e k (e ,l)- J Z(u;E )e k (E - u,l)du = . (5) 

o 

The boundary conditions for the moments are 

n(e ,Z)|, =0 = 0, 
?(£ ,0l*=o = 4 

If the relative energy loss of an electron in one collision is small the integro- 
differential equations (4) and (5) can be transformed by the Taylor expansion 
of integrands: 

d Id 2 
n(e - uj, I) ps n(e , 1) - uj — n(e , 1) + -uj 2 —^n(e , 1) , 



— — d — 1 d 2 — 

s k (e - uj, I) ps e k (e , 1) - u—e k (e , 1) + -uj 2 — 1 e k (e , 1). 



This gives the partial differential equations: 
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d did 

— n(e , 1) + (3( y e )—n( y e , 1) - 2 7 ( e °)fe2™( e °' = s ( £ o), (6) 



J^o, 1) + /3( £o ) A 5 k (5o? t) _ l 7(£o) JL 5 k (5o? /) = o. (7) 



The quantities £(eo), P( £ o), an d 7( £ o) in (6), (7) are the moments of the 
macroscopic differential cross-section: 



^max 

S(e ) = J T l (u;e )du, 



o 

^max 



/3(e ) = J ujY l (u;e )duj, 



o 

^max 



7(e ) = y cj 2 £(cj;e )du;. 
o 



The Eq. (7) for the second moment e 2 (e , 1) can be transformed into the equa- 
tion for the variance 

A(e ,l)=^(e ,l)-s 2 (e ,l). 

This equation is 



-A(e , 1) + P(e )—A{e , 1) - ^o)^ A{e , 1) 



7 (e )l— *(e ,J)l • (8) 



3 Cross-sections and related quantities 

The linear Compton scattering differential cross-section is 



dy x \ 1-y x(l-y) x 2 (l - y) 2 /' 
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e 2 

where y — — , a — nrL and r = - — — is the classical radius of electron. 

Eq [me 1 ) 1 

In the energy region of our interest the invariant dimensionless parameter x 
is small and the integral interaction coefficients S(e ), P(£o), and 7(^0) 

are described by the approximate formulas 



7 S 4 



p 



, _, 16 (mc 2 ) 2 

where L = — n L a , u; p = — . 

3 4a;o 

The quantities (5{sq) and 7(2:0) are the mean energy loss and the mean squared 
energy loss of an electron per unit path length. 



4 Solution of equations 



The partial differential equations (6), (7), and (8) can be solved by the iteration 
method. In the first approximation, where the terms with the second derivation 
are neglected, 



n 



[e , 1) = S / - 2 log + » E(e )J ( 1 + £„* (j^j J , (10) 



£0 



7 e 2 /3(e )l 
10a; p (l + P(e )l/e o y 



A(*o,0 = ^7? ,, , , (12) 



The second term in Eq. (10) is due to increasing of the interaction cross-section 
because of electron energy loss in collisions. But it is seen from the equation 
that this effect can be neglected in the energy region under consideration and 
the problem can be solved in one- velocity (X(e ) = const) approximation. 
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It follows from Eq. (12) that the variance A(e , has a maximum at the point, 
where 

P(eo)l = f (13) 

and 



e(go,0 _ 3 
e " 4 



Using Eqs. (12) and (11) one can derive the formula 

A( eo ,J) = ^-(l-- 

10 LU p \ Eq 



which was earlier obtained in [8]. 

It can be shown that in the one- velocity approximation the solution of Eq. (3) 
is the sum of the terms corresponding to nonscattered and scattered electrons: 

P(e|e , = exp (-S(e )Z) S(e - e) + P(e\e , 1) , 

oo 

P(e\e ,l) = Y,Pn(e ,l)U n (e\E ), 

71=1 

where 

(S(eo)O n 



Pn(£o,l) = exp(-E(e o )0 



is the Poisson distribution with the mean value n(e ,l) = S(e )/ and the 
function U n (e\e ) obeys the recurrent convolution formula 



w max 

U n (e\e ) = J w(u;eo)U n _i(e\e - u)du , 



where 

S(cj;e ) 



w(u;e ) 



S(e ) 



is the probability density function of u and 

Ui(e\e ) = w(e - e;e ) 
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It should be pointed out that U n (e\e ) is the energy distribution of n-scattered 
electrons. 

Decreasing of electrons energy due to their repeated collisions changes the 
energy spectrum of scattered photons. The photons produced in the secondary 
Compton scatterings are softer and the resulting spectrum can be written as 
the sum of terms corresponding to individual collisions: 



where $(u;|n) is the spectrum of photons resulting from the n— th electron 
collision and can be written in the form 



5 Numerical results 



The results of analytical calculations above agree with the data of our statis- 
tical simulation for 10 GeV electrons and 1 eV photon head-on collisions. It 
was supposed in this simulation that the number of electron collisions with 
laser photons is random. This number was selected from the Poisson distri- 
bution with fixed n. The simulation of individual collisions was carried out in 
the electron rest frame using the Klein and Nishina formula with the Lorentz 
transformation to the lab system. In the same way as in analytical calculations 
above we neglected the angular deflection of electrons but accounted for the 
energy decreasing after each collision. All results were obtained with statistics 
more than 10 6 trajectories. 

The energy spectra of electrons U n (e\eo) after fixed number of collisions are 
given in Fig. (1). The data on the electrons energy distributions P(e\e , 1) are 
given in Fig. (2) for several values of n. 



oo k 



k=l n=l 
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Fig. 1. Energy spectra of electrons after n = 1, 2, 3 collisions. Eq = WGeV, ojq = leV. 




It should be pointed out the discontinuous of the spectra for small n at the 
point £ = e Q — u max due to single scattered electrons and contribution of 
multiple scattered electrons at energies below the point of discontinuous. This 
contribution exists even for small n. 

It is seen from Fig. (2) that the width of the electron energy distributions 
decreases for such n which are greater than those one determined by Eq. (13). 

The energy spectra of photons <$>(u)\n) from the n— th scattering of electron 
and the resulting spectra for several n are shown in Fig. (3) and Fig. (4). It is 
seen from Fig. (4) that the contribution of multiple scattering should be taken 
into account even for the photon target of small thickness. 
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Fig. 3. Energy spectra of photons generated in the n-th collision, n = 1, 2, 3. 
e = lOGeV, oj = leV. 
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Fig. 4. Energy spectra of photons for n = 0.5, 1. Eq = lOGeV, u>o = leV. 
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Fig. 5. Energy spectra of photons with emission angle (9 7 < — — — for n = 0.5, 1, 



2. e = WGeV, uj = leV. 



2 £ 



Fig. (5) shows that the electrons multiple collisions influence on the spectra 
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1 TfleP' 

of photons with emission angle 9 1 < . It is seen from the figure that 

2 Eq 

multiple Compton scattering results in broadening of the photon spectra with 
increasing of the photon target thickness. 

The mean energy loss of electrons in backward Compton scattering and the 
variance of the electron energy distribution are shown in Fig. (6) and Fig. (7). 
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Fig. 6. Mean energy of electron in back Compton scattering, e$ = WGeV, ujq = leV. 
Points - simulation, solid line - Eq. (11). 
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Fig. 7. Variance of electron energy distribution, eo = WGeV, ujq = leV. Points 
simulation, solid line - Eq. (12). 



It is seen that the Monte Carlo data agree with the analytical calculations 
above. 
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